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Progress by the Lattice Hadron Physics Collaboration in determining the baryon 
and meson resonance spectrum of QCD using Monte Carlo methods with space- 
time lattices is described. The extraction of excited-state energies necessitates 
the evaluation of correlation matrices of sets of operators, and the importance of 
extended three-quark operators to capture both the radial and orbital structures 
of baryons is emphasized. The use of both quark-field smearing and link-field 
smearing in the operators is essential for reducing the couplings of the operators 
to the high-frequency modes and for reducing statistical noise in the correlators. 
The extraction of nine energy levels in a given symmetry channel is demonstrated, 
and identifying the continuum spin quantum numbers of the levels is discussed. 



1. Introduction 

A charge from the late Nathan Isgur to use Monte Carlo methods to ex- 
tract the spectrum of hadron resonances and hadronic properties resulted in 
the formation of the Lattice Hadron Physics Collaboration (LHPC) in the 
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year 2000. As part of a national collaboration of lattice QCD theorists, the 
LHPC acquired funding from several sources, including the DOE's Scientific 
Discovery through Advanced Computing initiative, to build large comput- 
ing clusters at JLab, Fermilab, and Brookhaven, as well as to develop the 
software to carry out the needed large-scale computations. The LHPC has 
several broad goals: to compute the spectrum of QCD from first principles, 
to investigate hadron structure by computing form factors, structure func- 
tions, and other matrix elements, and to study hadron-hadron interactions. 
This talk focuses solely on our efforts to determine the spectrum of QCD. 

Extracting the spectrum of resonances in QCD is a big challenge. The 
determination of excited-state energies requires the use of correlation ma- 
trices of sets of operators. The masses and widths of unstable hadrons 
(resonances) must be deduced from the energies of single-particle and multi- 
hadron stationary states in a finite-sized box. This necessitates the use of 
multi-hadron operators in the correlation matrices, and the computations 
must be performed in full QCD at realistically-light quark masses. For 
these reasons, the computation of the QCD spectrum is indeed a long-term 
project. This talk is a brief status report of our efforts. 

I first describe in Sec. 2 how excited-state energies are extracted in our 
Monte Carlo calculations, and discuss issues related to unstable resonances. 
Our construction of baryon operators is then outlined in Sec. 3 (see Ref. 1 
for further details), and in Sec. 4, the importance of using smeared quark 
and gluon fields is demonstrated, (see also Ref. 2). The first-time extraction 
of nine nucleon energy levels is also presented in this section. Concluding 
remarks are made and future work is outlined in the concluding Sec. 5. 

2. Excited states and resonances 

In the path integral formulation of quantum field theory with imaginary 
time, stationary-state energies are extracted from the asymptotic decay 
rates of temporal correlations of the field operators. If <I>(i) is a Hcisenberg- 
picture operator which annihilates the hadron of interest at time t, then its 
evolution <j>(t) = e Ht Q(0)e~ Ht , where H is the Hamiltonian, can be used 
to show that in a finite-sized box, 



C(t) = <0|$(t)$t(0)|0) = ^(0|e m $(0)e- fft |n)(n|$t(0)|0), (1) 
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where {\n}} is the complete set of discrete eigenvectors of H. We assume 
the existence of a transfer matrix, and temporal boundary conditions have 
been ignored for illustrative purposes. One can then extract A\ and E\ — Eo 
as t ^ oo, assuming (0|$(0)|0) = and (0|$(0)|1) ^ 0. 

A convenient visual tool for demonstrating energy-level extraction is the 
so-called "effective mass" defined by m c s{t) — ln[C(t)/C(t + a t )], where t is 
time and a t is the temporal lattice spacing. The effective mass tends to the 
actual mass (or energy) of the ground state as t becomes large, signalled 
by a plateau in the effective mass. At smaller times before this plateau is 
observed, the effective mass varies due to contributions from other states 
in the spectrum. The correlation C{t) is estimated with some statistical 
uncertainty since the Monte Carlo method is used, and usually, the ratio 
of the noise to the signal increases with temporal separation t. Hence, 
judiciously-chosen operators having reduced couplings with contaminating 
higher-lying states are important in order to observe a plateau in the effec- 
tive mass before noise swamps the signal. Key ingredients in constructing 
such operators are the use of smeared quark and gluon fields, the incorpo- 
ration of spatially-extended assemblages of the fields, and the use of sets of 
different operators to exploit improvements from variational methods. 

Methods of extracting excited-state energies are well known 3 ' 4 . For a 
given N x N Hermitian matrix of correlations C a p(t) = (0|$ ct (t)$^(0)|0}, 
the N principal correlators X a (t,t n ) are defined as the eigenvalues of the 
matrix C(to)~ 1 / 2 C{t)C(to)~ 1 / 2 , where t a is some reference time (typically 
small), and one can show that 

lim A a (t,to) = e-(*- t0 ^(f+O(e-* A ^)), AE a = min \E P -E a \, (3) 

assuming Eq = and Ai > A2 > A3 • • • . The N principal effective masses 
mjf (() = ln[A a (t, t n )/X a (t+a t , to)} now tend (plateau) to the N lowest-lying 
stationary-state energies which couple to the N operators. The associated 
eigenvectors are orthogonal, and a knowledge of the eigenvectors can yield 
information about the partonic structure of the states. 

Many of the hadron states we wish to study are unstable resonances. 
Our computations are done out of necessity in a box of finite volume with 
periodic boundary conditions. Hence, the momenta of the particles we 
study are quantized, so all states are discrete in our computations. Thus, 
we can only determine the discrete energy spectrum of stationary states in 
a periodic box, which are admixtures of single hadrons and multi-hadron 
states. Resonance masses and widths must somehow be deduced from the 
finite-box spectrum 5 ' 6 ' 7 ' 8 . Once the masses of the stable single particle 
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Table 1. Continuum limit spin identification: the number n^ of 
times that the A irrep. of the octahedral point group Oh occurs in 
the (reducible) subduction of the J irrep. of SU{2). The numbers for 
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states have been determined, the placement and pattern of their scattering 
states are known approximately, and the dependences of their energies on 
the volume are roughly known. Resonances show up as extra states with 
little volume dependence. Our initial goal is simply to ferret out these 
resonances, not to pin down their properties to high precision. 

3. Operator construction 

Our approach to constructing hadron operators is to directly combine the 
physical characteristics of hadrons with the symmetries of the lattice reg- 
ularization of QCD used in simulations. For baryons at rest, our operators 
are formed using group-theoretical projections onto the irreducible repre- 
sentations (irreps) of the Oh symmetry group of a three-dimensional cubic 
lattice. There are four two-dimensional irreps G\ g , G\ u , G?2 g , Gi u and two 
four-dimensional representations H g and H u . The continuum-limit spins J 
of our states must be deduced by examining degeneracy patterns across the 
different Oh irreps (see Table 1). For example, a J p = j + state will show 
up in the G\ g channel without degenerate partners in the other channels, 
and a J p = | state will show up in the H g channel without degenerate 
partners in the other channels. Four of the six polarizations of a J p = | + 
state show up as a level in the H g channel, and the other two will occur as a 
degenerate partner in the G2 g channel, whereas three degenerate levels, one 
in each of the three G\ g , G2 g , H g channels, may indicate a single J p = | + 
state (or the accidental degeneracy of a spin-i and a | state). 

Baryons are expected to be rather large objects, so local operators will 
not suffice. Our approach to constructing spatially-extended operators is to 
use covariant displacements of the quark fields along the links of the lattice. 
Displacements in different directions are used to build up the appropriate 
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single- singly- doubly- doubly- triply- triply- 
site displaced displaced-I displaced-L displaced- T displaced-0 

Figure 1. The spatial arrangments of the extended three-quark baryon operators used. 
Smeared quark fields arc shown by solid circles, line segments indicate gauge-covariant 
displacements, and each hollow circle indicates the location of a Levi-Civita color cou- 
pling. For simplicity, all displacements have the same length in an operator. 

orbital structure, and displacements of different lengths can build up the 
needed radial structure. All our three-quark baryon operators are superpo- 
sitions of gauge- invariant, translationally- invariant terms of the form 

^ 7 %(t) =5>fc= (D^{x,t))t (D^(x,t))^ {D^(x,t))% (4) 

where A,B,C indicate quark flavor, a, b, c are color indices, a, ft, 7 are 
Dirac spin indices, tp indicates a smeared quark field, and denotes the 
smeared p-link covariant displacement operator in the j-th direction. The 
smearing of the quark and gauge field will be discussed later. There are 
six different spatial orientations that we use, shown in Fig 1. The singly- 
displaced operators are meant to mock up a diquark-quark coupling, and 
the doubly-displaced and triply-displaced operators are chosen since they 
favor the A-flux and F-flux configurations, respectively. 

Next, the ^afi-^ijk are combined into elemental operators B^(t) having 
the appropriate flavor structure characterized by isospin, strangeness, etc. 
We work in the m u = (equal u and d quark masses) approximation, 
and thus, require that the elemental operators have definite isospin, that is, 
they satisfy appropriate commutation relations with the isospin operators 
T3, r+, t_. Since we plan to compute full correlation matrices, we need not 
be concerned with forming operators according to an ST/ (3) flavor symme- 
try Maple code which manipulates Grassmann fields was used to identify 
maximal sets of linearly independent elemental operators. The operators 
used are shown in Table 2, and the numbers of such independent operators 
are listed in this table as well. 

The final step in our operator construction is to apply group-theoretical 
projections to obtain operators which transform irreducibly under all lattice 
rotation and reflection symmetries: 

B AXF {t)= ifA_ £ r (A) (jR) Ur B r {t) jji (5) 
RtZ h 
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Table 2. Elemental operators for various baryons (left). The numbers of linearly inde- 
pendent elemental operators of each spatial kind (right). 



Baryon 



Operator 



Spatial Type 



N 



A 



A++ 
S+ 
N+ 
5° 
A° 



ocfi~f,ijk 

QSSU 

a(3^,ijk 
a(3'y,ijk ap^,ijk 

QSSS 

ap~/,ijk 



single-site 20 20 40 24 

singly-displaced 240 384 624 528 

doubly-displaced-I 192 384 576 576 

doubly-displaccd-L 768 1536 2304 2304 

triply-displaccd-T 768 1536 2304 2304 

triply-displaced-O 512 1024 1536 1536 



where A refers to an Oh irrcp, A is the irrep row, go h is the number of 
elements in Oh, g?a is the dimension of the A irrep, Tmn(R) is a A represen- 
tation matrix corresponding to group element R, and Ur is the quantum 
operator which implements the symmetry operations. The projections in 
Eq. (5) are carried out using computer software written in the Maple 9 sym- 
bolic manipulation language. 



4. Field smearing and operator pruning 

For single-site (local) hadron operators, it is well known that the use of 
spatially-smeared quark fields is crucial. For extended baryon operators, 
one expects quark-field smearing to be equally important, but the relevance 
and interplay of link-field smearing is less well known. Thus, we decided 
that a systematic study of both quark-field and link- variable smearing was 
warranted. 

The link variables were smeared U — > U using the analytic stout link 
method of Rcf. 10. There are two tunable parameters, the number of it- 
erations n p and the staple weight p. For the quark-field, we employed 
gauge-covariant Gaussian smearing 

*(z)=(l + ^A)%(z), (6) 
A*(a;)= J2 (Uk(xMx + k)-*{xj), (7) 

fe=±l,±2,±3 

where A denotes the smeared three-dimensional gauge-covariant Laplacian. 
The two parameters to tune in this smearing procedure are the smearing 
radius a s and the integer number of iterations n a . 
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Figure 2. Effective masses M(t) for unsmeared (circles) and smeared (triangles) opera- 
tors Ogs, Osd, Otdt- Top row: only quark- field smearing n CT = 32, a s = 4.0 is used. 
Middle row: only link-variable smearing n p = 16, n p p = 2.5 is applied. Bottom row: 
both quark and link smearing n CT = 32, a s = 4.0, n p = 16, n p p = 2.5 are used, dra- 
matically improving the signal for all three operators. Results are based on 50 quenched 
configurations on a 12 3 X 48 anisotropic lattice using the Wilson action with a s ~ 0.1 
fm, a s /at ~ 3.0. The quark mass was chosen so that ra^ ~ 700 MeV. 



For our tests of the efficacy of quark-field and gauge-link smearing, 
correlators were computed for three particular nucleon operators: a single- 
site operator Oss in the G\ g irreducible representation of the cubic point 
group, a singly-displaced operator Osd with a particular choice of each 
Dirac index, and a triply-displaced-T operator Otdt with a specific choice 
of each Dirac index. Our findings are summarized in Fig. 2. The top row 
shows that applying only quark-field smearing to the three selected nucleon 
operators significantly reduces couplings to higher-lying states, but the dis- 



February 1, 2008 5:21 Proceedings Trim Size: 9in x 6in 



Nstar2005'writeup 



1 

^0.8 
"0.6 
«** 0.4 

0.2 




"0.6 
«** 0.4 

0.2 


1 

^0.8 
"0.6 
«** 0.4 
0.2 




I 1 1 1 1 I 1 1 1 

ground state 



' 1st' excited state' 



2nd' excited' state' 



■ ■ I ■ I ■ 



i i i I i i i i I i i i i 

3rd excited state . 



■ 1 1 1 1 1 1 1 1 1 1 1 1 



i i i I i i i i I i i i i 

4th excited state 



■ 1 1 1 1 1 1 1 1 1 1 1 1 



I I | I I I I | I I I I 

5th excited state . 



"••••••hj.i;" 



■ ' 1 ■ ' ' ' 1 ' ' ' ' 



1 1 1 1 1 1 1 1 1 1 1 1 1 

6th excited state . 



•••■55 



J 1 



. .. i .. .. i . .. . 



. 1 1 1 1 1 1 1 1 1 1 1 1 1 . 

■_ 7th excited state j 






■ 1 1 1 1 1 1 1 1 1 


III' 



I I I I I I I I I I I I 

8th excited state . 



7 







5 10 

tla. 



15 



5 10 

t/a 



15 



10 



t/a 



15 



Figure 3. Principal effective masses for the lowest-lying nine states in the Gi g nuclcon 
channel. Results are based on 100 quenched configurations on a 12 3 X 48 anisotropic 
lattice using the Wilson action with a s ~ 0.1 fm, a s jat ~ 3.0. The quark mass was 
chosen so that m n ~ 700 MeV. Quark-field smearing n a = 32, a s = 3.0 and link- 
smearing rip = 16, rtpp = 2.5 are used. These principal effective masses come from 
a 10 X 10 correlation matrix, where 2 singly-displaced, 5 doubly-displaccd-I, 2 doubly- 
displaccd-L, and 1 triply-displaccd-T operators were chosen. 



placed operators remain excessively noisy. The second row illustrates that 
including only link-field smearing substantially reduces the noise, but does 
not appreciably alter the effective masses themselves. The bottom row 
shows dramatic improvement from reduced couplings to excited states and 
dramatically reduced noise when both quark-field and link-field smearing 
is applied, especially for the extended operators. The "best" quark-field 
smearing parameters n a and a were determined by requiring that the ef- 
fective mass for the three chosen operators reach a plateau as close to the 
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source as possible. The gauge-link smearing parameters were tuned so as 
to minimize the noise in the effective masses. One interesting point we 
also learned was that the preferred link-smearing parameters determined 
from the static quark-antiquark potential produced the smallest error in 
the extended baryon operators. 

The computation of correlation matrices using hundreds of operators is 
not feasible, so it is necessary to "prune" out unnecessary operators. The 
first step in this pruning is to examine the effective masses of the diagonal 
elements of the correlation matrices to identify and eliminate noisy oper- 
ators. Keeping only operators with small statistical uncertainties yields a 
set of about forty to fifty operators in each symmetry channel. We then 
computed the correlation matrix of this reduced set of operators, examin- 
ing whether further reductions to the operator sets could be made without 
increased contamination in the principal effective masses. These computa- 
tions are still ongoing, but preliminary results are shown in Fig. 3. This 
figure shows that it is possible to extract at least nine levels in a given sym- 
metry channel, a feat which has never before been accomplished. Demon- 
strating that this number of energy levels can be reliably extracted is an 
important milestone in our long-term project. 

5. Conclusion 

We have outlined a program to study the resonance spectrum in lattice 
QCD. The use of the variational method and the need to isolate several 
energy levels in each channel require a sufficiently broad basis of operators. 
Having developed suitable group-theory methods to project operators onto 
the irreducible representations of the cubic group, and having examined 
the efficacy of both quark- and gaugc-link-smearing, we are now identify- 
ing a more limited set of operators that we will employ in a large-scale 
study of the hadron spectrum. Our methods are applicable not only to 
baryons, but also to mesons, tetra-quark and pentaquark systems, and to 
states with excited glue. Only by performing such a program can we hope 
to identify the states of QCD, and in particular their spins and parities, in 
the continuum limit. Ultimately, when quark loops are included at real- 
istically light quark masses, multi-hadron (baryon-meson) operators must 
be included in our correlation matrices, and finite-volume techniques will 
need to be employed to ferret out the baryon resonances from uninteresting 
scattering states. We are currently exploring different ways of building such 
operators. 
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